Feshbach-resonance-induced atomic filamentation and quantum pair correlation in 

atom-laser-beam propagation 
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We study the propagation of an atom laser beam through a spatial region with a magnetic field 
tuned to a Feshbach resonance. Tuning the magnetic field below the resonance produces an effective 
focusing Kerr medium that causes a modulational instability of the atomic beam. Under appropriate 
circumstances, this results in beam breakup and filamentation seeded by quasi-particle fluctuations, 
and in the generation of correlated atomic pairs. 
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It has been recognized for quite some time that 
the two-body interactions that govern the dynamics of 
atomic Bose-Einstein condensates (BECs) , and which are 
expressed most simply in terms of the s-wave scattering 
length a, play a role for matter waves analogous to that 
of a cubic (Kerr) nonlinearity in optics. This was illus- 
trated in a number of experiments that include matter- 
wave four- wave mixing |jj , the generation of dark 0] and 
bright atomic solitons, coherent matter- wave ampli- 
fication and the generation of quantum-correlated 
atomic pairs H. It has also been demonstrated that the 
strength of this nonlinearity can be varied, and its sign 
changed, using in particular Feshbach resonances || 0. 
In the language of nonlinear optics, changing the sign of 
the scattering length amounts to switching from a defo- 
cusing (for repulsive interactions, a > 0) to a focusing 
(for attractive interactions, a < 0) nonlinear medium. 
The tunability of the scattering length offers a degree 
of flexibility of considerable interest in the study of the 
manybody properties of atomic condensates, as well as 
for nonlinear atom optics applications. 

As pointed out by Timmermans et al. S, caution is 
in order when considering the dynamics of a condensate 
in the vicinity of a Feshbach resonance. In the atomic 
systems under consideration, they correspond to a situa- 
tion where the total energy of two colliding atoms equals 
the energy of a bound molecular state. In the vicinity 
of the resonance, it is essential to properly account for 
the molecular dynamics. Indeed, it has recently been 
unambiguously established experimentally that a proper 
description of the system requires that one accounts for 
the coherent coupling between the atomic and molecular 
field, resulting in the onset of atom-molecule coherence 
in BECs ||. The presence of the molecular state is also 
predicted to lead to "resonant superfluidity" in quantum- 
degenerate atomic Fermi systems [ flO| . 

Far from resonance, however, it is possible to adia- 
batically eliminate the molecular field from the system's 
dynamics, in which case the atomic condensate is subject 
to an effective, tunable scattering length that describes 
the binary atomic collisions. The goal of this letter is to 
illustrate how this tunable scattering length can be ex- 



ploited to manipulate and control the propagation of a 
continuous-wave atom laser beam. The specific geom- 
etry that we consider consists of an atom laser beam 
propagating through a region where a spatially depen- 
dent magnetic field is applied to switch the nature of the 
binary collisions from being repulsive, or defocusing, to 
attractive, or focusing. This results in a spatial change 
of the interatomic interaction similar to the presence of 
a Kerr medium in nonlinear optics, albeit with some im- 
portant differences. We show that under appropriate cir- 
cumstances, the quasi-particle fluctuations in the initial 
atom laser beam can seed the filamentation of the atom 
laser beam, with the generation of pair-correlated atoms 
in spatial sidebands. This is similar to the amplification 
of squeezed light by parametric processes in a nonlinear 
optical Kerr medium, but with the important difference 
that in the present situation it is not necessary to start 
from a squeezed input. 

We recall that the generation of correlated atomic pairs 
from BECs has recently been studied by a number of au- 
thors ||, O]. In these works, correlated atomic pairs are 
generated cither from vacuum fluctuations [pd| , or grow 
via four-wave mixing from the combined effects of a seed 
beam and vacuum fluctuations Js). In contrast, the pro- 
posed scheme involves the generation of a pair-correlated 
atomic source from quasi-particle vacuum fluctuations, 
with an atom laser beam as the sole pump beam. 

A Feshbach resonance is a scattering resonance where 
two particles in an open channel collide to form a bound 
state in the closed channel. In the case of quantum- 
degenerate atomic fields, such as Bose-Einstein conden- 
sates or atom laser beams, the bound state in the closed 
channel is a coherent molecular state. The dynamical 
equations describing the atom-molecule coherent cou- 
pling can be written as [B[ 
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Here $ a . m are the atomic and molecular field opera- 
tors, m is the atomic mass, U a = (Airh 2 /m)a a , U m = 
(2ttTl 2 /m)a m , and U am — (3nh 2 j1m)a am are the atom- 
atom, molecule-molecule, and atom-molecule binary in- 
teraction strengths, with a a , a m , and a am the cor- 
responding scattering lengths, and a is the coupling 
strength the Feshbach resonance. The parameter e is 
the energy difference between the molecular state and 
the atomic state. 

For large resonance detunings |e|, it is possible to adi- 
abatically eliminate the molecular field, leading to the 
effective Heisenberg equation of motion 

+ **!*•)*•■ < 2 > 

for the atomic Schrodinger field. In this expression, the 
effective coupling constant h\ = U a — 2a 2 /e accounts 
for the change of the the strength of the atomic binary 
collisions resulting from the Feshbach resonance. 

We now proceed to study the propagation of an atom 
laser through a spatially-dependent effective nonlinearity, 
which for simplicity we take to be of the form 

h X (z) = U a -2(a 2 /e)0{z), 

where &(z) is a Heaviside step function. In particular, 
we consider the situation in which the binary collisions 
change from repulsive for z < 0, to attractive for z > 0. 
The continuous wave atom laser beam is incident from 
z < 0, and we describe the incident atomic field as a co- 
herent beam of Bose-condensed atoms propagating along 
the z-axis with wave vector fcoz = (mva/K)i. We also 
ignore any atoms that are back reflected at the interface: 
If we take account of the fact that the magnetic field ac- 
tually varies over a length L around z — 0, then using 
a WKB analysis we find that back reflections are negli- 
gible if (27r/fco)(2a 2 /|e|)no/£ <C v 2 /2m where no is the 
density of atoms in the condensate. 

An important distinction between this and the generic 
situation considered in nonlinear optics, where the beam 
is incident from vacuum, is that the incident atomic field 
is self-interacting. Hence, a proper description must in- 
clude the quasi-particle vacuum fluctuations propagating 
together with the coherent BEC component. As such, a 
realistic atom laser beam consists of the sum of a moving 
condensed component, described classically, and quan- 
tum fluctuations ]i2| , 

$ a (r, t) = [yfa+ </>(r, f)] e «*o-*»)*-*"o*, z < 0, 

_ (3) 
where y/no describes the coherent BEC component 
(taken real without loss of generality), tuvo = h 2 k 2 ./2m, 
and k n = x n o/vo is a nonlinear phase shift. The noise 
operator ip(r,t) obeys bosonic commutation relations. 



From the corresponding Heisenberg equation of motion 
, one observes that the binary collisions contribute a 
term of the form 

lh ^W^ ~ UaU ° (^ (r ' t] + ^ t(r ' t] ) 

to the evolution of the noise operator ip(r,t). This form of 
evolution is quite familiar in quantum optics, where it is 
known to lead to the generation of squeezed states. Thus, 
the repulsive self-interaction of the incident atom laser 
beam for z < is sufficient to provide initial squeezed 
state vacuum fluctuations at the input port of the (fo- 
cusing) nonlinear medium z = 0. 

We proceed by expanding the noise operator "0(r, t) in 
terms of plane-waves as 

$(r,t) = JdqJ dk z [a (<i,k z ,ty k ^ z - v ^ 

-6o(qA,*)e- <Mz ^ o4) ]e icl - rx , (4) 

where q is the transverse wave vector and the operators 
fio(q, k z ,t) and bo(q 7 k z ,t) = — ao(q> — k z ,t) describe for- 
ward and backward propagating modes in the rest frame 
of the atomic beam, respectively. 

In the incident region with repulsive interactions, z < 
0, it is convenient to express the small quantum fluctu- 
ations about the condensate wave function in terms of 
quasi-particles characterized by bosonic annihilation op- 
erators d(q, fc z , t) and /3(q, k z ,t) = — a(— q, —k Z} t), given 

by @ 

fio(q, k z ,t) = Wfea(q, k z ,t) +Wfc/3 t (q, k z ,i), 
h (q,k z ,t) = u fe /3(-q, k z ,t) +v fc d t (-q, k z ,t), 

where 
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with k = ^q 2 + k z , et- = h 2 k 2 /2m, and Ek — 
yj e'f. + 2ek^X n o being the quasiparticle energy associated 
with d(q, k x , t) and /3(q, k z ,t). 

We are now in a position to study the propagation of 
the atom laser beam as it enters the region of Feshbach 
resonance, z > 0. Specifically, in this paper we focus 
on the initiation of the filamentation process before ap- 
preciable depletion. For this purpose, it is sufficient to 
consider a linearized perturbation theory in the slowly 
varying amplitude approximation. For z > we decom- 
pose the atomic field operator in the same manner as 
Eq. (|J) with the substitution of ^(r, t) in place of ^(r, t) 
for the fluctuation operator. Furthermore, analogous to 
Eq. (U), we also decompose the noise operator $(r, t) in 
terms of quasi-plane waves, 

$(r,t) = / dq$ q (z,i)e iqr -, (5) 
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Here we allow for a slow z-dependence of the operators 
& and b. This approach, familiar from quantum optics, 



allows one to account for the possible amplification of the 
quantum noise. Spatial variation of the quasi-particle 
operators was ignored in the region z < due to the 
absence of sideband gain for repulsive binary collisions. 

Substituting Eq. (||) into Eq. (||) yields pairs of coupled 
operator equations that describe the linearized spatial 
evolution of the condensate sidebands, 
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where T = t + z/vq. Assuming a solution to Eq. (Q) of the form e gT and matching the field operator solutions in the 
two regions at z = we obtain finally 



a(q,k z ,z) = {[U(q,z)u k + V(q, z)v k ] a(q, k z ) + [U(q,z)v k + V(q, z)u k ] /3 f (q, k z )\, 
S(q, k z , z) = - j [V(q, z)u k + U (q, z)v k ] a f (-q, k z ) + [V(q, z)v k + U (q, z)u k ] /3(-q, k z ) } , 
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where the coefficients U and V are given by 



U(q, z) — cosh I g 
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V(q,z) 
with the gain parameter 



9 = V (X"o) 2 - (e 9 + X"o) 



(9) 



(10) 



Without loss of generality we have taken t = in order to 
focus on how the filamentation develops with propagation 
distance. 

We observe that the linearized solution of Eqs. (||) can 
exhibit gain in the case of a focusing Feshbach resonance, 
X < 0. Specifically, the gain coefficient g is real for q < 
1\J m\x\no/ft, with a maximum, g max — \x\ n o, for the 
transverse wavevector 



y/2m\x\n /h. 



(11) 



Physically, this growth of a range of coupled pairs of spa- 
tial sidebands ±q around the condensate corresponds to 
a transverse modulational instability with a cone of un- 
stable transverse momenta peaked around |q| = <7 max - 
If allowed to grow this instability will ultimately lead 
to the breakup of the incident laser beam into trans- 
verse structures or filaments similarly to the situation 
in nonlinear optics |l^, If instead of using a plane- 
wave condensate we consider a beam with widths w Xj y 
in the x and y transverse directions, we can eliminate 
the instability in the y-direction, for example, by choos- 



wavevector cannot be supported in the y-direction. In 
that case the cone of unstable wavevectors becomes a 
pair of sidebands that are peaked at ±<? max along the x- 
direction. This quasi-one-dimcnsional geometry may be 
better to isolate correlated atomic pairs. 

This modulational instability is intimately related to 
the collapse of 85 Rb condensates seen in recent experi- 
ments 0. In these experiments, stable repulsive BEC's 
with densities of no w 10 12 cm~ 3 are produced in the 
presence of a magnetic field tuned close to a Feshbach 
resonance (B « 160 G) and then the magnetic field is 
tuned above 167 G where the interaction becomes attrac- 
tive. Far from the resonance the interactions are attrac- 
tive with a background scattering length of a^g = —23.85 
nm jfj). Using these parameters gives a <7 max = 7.75 x 10 6 
m _1 or a transverse velocity of 0.6 cm/s for the sideband 
atoms. 

Proceeding with the example of a plane-wave conden- 
sate, to quantify the statistical properties of the gen- 
erated atomic sidebands, we evaluate the field correla- 
tion and intensity correlation functions for &q(z) defined 
in Eq. (^). Assuming that the quantum fluctuations in 
the incident atom laser beam are in the quasi-particle 
vacuum state l^o)) i- e -> Q^o) = $\^o) — 0; an d with 
Eqs. (|) and (§) we find, 



(m^M = p(q)p(q') + \Q(Q)\ 2 Sc l ,-c l ', 



(12) 
(13) 



where all the operators are evaluated at the same longi- 



tudinal distance z and in Eq. (13), we have taken q^q' 



ing q max w y < l,q n 



3> 1, so that the most unstable Physically, p(q) is the probability density in momentum 



4 



space of finding an atom with a transverse momentum 
?iq after propagating a distance z in the Feshbach re- 
gion, and \Q(q)\ 2 is the joint probability that if one atom 



is detected with momentum tiq a second one will be de- 
tected with opposite momentum —Tic\. They are given 
by, 



Q(q) 
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FIG. 1: Correlation spectrum of the correlated atomic pairs 
calculated from Eq. (|l4|) . Due to the divergence at large wave 
number, the upper limit of k z is cut off at k s = mcp/h — the 
wave number corresponding to the Bogoliubov sound velocity. 
Here the wave numbers, energy and length are in units of k a , 
h 2 k 2 1 (2m) and l/k s , respectively. The parameters used in the 
calculation are vo = 5c p , x n o = —0.5. The plots show that 
both p and \Q\ grow fastest near g max = I/a/3 [see Eq. (Ell)]. 

with e q = e q + xn . 

As an illustration we assume that atoms are detected 
after having propagated a distance z in the Feshbach 
region such that the linearized analysis is still valid. 
The rate at which atoms are detected at an angle 9 — 
tan (q/ko) relative to the z-axis will be proportional to 
p(q) while the rate of coincidence detections at +9 and —9 
will be proportional to |Q(q)| 2 . The corresponding prob- 
ability densities are plotted in Fig. [j], from which we see 
that, as expected physically, with increasing propagation 
distance z the most probable transverse wavevector for 
detecting either an off-axis atom or a correlated atomic 



pair at ±q is increasingly |q| = q max . These results show 
that the transverse modulational instability can indeed 
be an effective source of correlated atomic pairs. 

In conclusion, we have studied the propagation of an 
atom laser beam through a Feshbach resonance region 
where the effective interaction between atoms is tuned 
from being repulsive to being attractive. The attractive 
interaction is equivalent to a self-focusing Kerr medium 
for the atom laser beam which can generate quantum- 
pair-correlated atomic filament beams by a modulational 
instability. This opens up new window to create non- 
classical ultracold atomic sources by coherent control of 
the atomic interactions. The quantum-pair-correlated 
atomic beam may have potential applications in both 
fundamental research and applications in quantum in- 
formation and computation. 
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